ESTIMATION OF A COST FUNCTION FOR A 
NAVAL AIR REWORK FACILITY 



Wilbur Cobb Trafton 



Library 

Naval Postgraduate School 
Monterey. California 93940 



U iff H P 

Mill. 




nterey, bam 



5 p p p ft 

UfiUU 




H l^sss^l ^jbsj v&sas?^ Sj ^--r^ 



ESTIMATION OF A COST FUNCTION FOR A 
NAVAL AIR REWORK FACILITY 



by 



Wilbur Cobb Trafton 



Thesis Advisor: 



Norman K. Womer 



March 1973 



T 154782 



App/Lovzd j fiosi pubtic Acte.ai e; cUzt'uJnition unZimLtzd . 



Estimation of a Cost Function for a 
Naval Air Rework Facility 



by 



Wilbur Cobb Trafton 
Lieutenant, United States Navy 
B, S. , United States Naval Academy, 1966 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
March 1973 






' 



/ 




Library 

Naval Postgraduate School 
Monterey, California 93940 



ABSTRACT 

The objective of this study was to estimate a cost function from 
a Constant Elasticity of Substitution production function and a Cobb- 
Douglas production function for the aircraft rework and engine repair 
programs at the Naval Air Rework Facility, North Island, San Diego, 
California, The cost functions were estimated by multiple regression 
analysis , from data aggregated from actual data taken from production 
records of the two programs. An attempt was made to validate the two 
cost functions that were obtained, and a methodology was outlined for 
comparing predicted costs to actual production costs at the Naval Air 
Rework Facility, 
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I. INTRODUCTION 



A. ORGANIZATION, NARFNI 

The Naval Air Rework Facility, Naval Air Station, North Island, 

San Diego, California, (NARFNI) is one of seven rework facilities serving 
the U„S. Navy and the U.S. Marine Corps. It is directly responsible for 
all major maintenance, incorporation of technical changes, and repair of 
extensively damaged West Coast based F-4, F~8, C-2, and E-2 aircraft and 
H-3, H-46, and H-53 helicopters. 

Bradley [l] covers the organization and operating proceedures at 
NARFNI in a very detailed manner. What follows is a brief description 
of NARFNI. 

NARFNI is directly responsible to the Naval Air Systems Command 
Representative, Pacific, who through the Commander, Naval Air Systems 
Command, is responsible to provide . aircraft maintenance to the Commander- 
in-chief, Pacific Fleet, in the operational chain of command. 

NARFNI is composed of nine major divisions (Figure 1), and each 
division is composed of a direct labor force and an indirect labor force. 
The direct labor force encompasses all skilled tradesmen, while the 
indirect labor force encompasses all managerial and administrative per- 
sonnel , 

All aircraft in the Naval Service are scheduled for preventive 
maintenance on a regular Progressive Aircraft Rework Cycle (PAR Cycle) 
during their life span. While undergoing PAR, necessary modifications 
to existing systems are made and major overhaul of component parts is 
carried out. Prior to induction into a specific PAR Cycle, an estimate, 
based on historical data of man hours required to rework similar aircraft, 
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is obtained. This estimate is referred to as Production Load Norm (NORM). 
Managers then plan the rework cycle for each aircraft, ordering material 
as required and allocating existing manpower, 

NARFNI is engaged in three major programs: 

1. Airframe Repair 

2. Engine Repair 

3. Component .Repair 

Most of the work is for the Naval Air Systems Command Representative, 
Pacific. Small amounts of work are done for other customers, particularly 
the Air Force and the Coast Guard. 

Most of the work at NARFNI is done on a fixed price contract basis. 
Costs negotiated in this contract between the customer and NARFNI include 
direct labor and" direct material costs, as well as an overhead cost to 
cover indirect labor, station support services, and employee benefits. 

NARFNI is funded primarily by the Navy Industrial Fund (NIF), 
Basically the NIF system provides a working capital fund to finance 
continuing cycles of ongoing operations. Receipts derived from these 
operations are used to perpetuate the fund. All work done by NARFNI is 
paid for initially out of the NIF which is then reimbursed from customer 
appropriations. NARFNI is required by NAVAIRSYSCOM to plan and control 
its financial affairs in an attempt to incur zero loss at the end of 
each fiscal year. Prices for services rendered during each quarter are 
established through negotiations with the customer at the beginning of 
each quarter during the Quarterly Fleet Readiness Support Conference. It 
is NARFNI 1 s goal to balance deficits or surplusscs in one quarter with 
compensating surplusses and deficits the next, thus maintaining the NIF 
at a constant value. 



9 



B. WIPICS 



In 1972 NARFNI installed a computerized Work in Proeess Inventory 
Control . System (WIPICS) developed by the Rohr Corporation, Chula Vista, 
California, Spooner [2] gives a good overview of the history behind the 
installation of WIPICS at NARFNI and the details involved in actual day 
to day operation of WIPICS. The system employs many of the latest advances 
in eomputer technology which include real time data file updating and 
information retrieval with communications via Touch Tone^(R ; telephones 
and remote teletypewri te terminals. Basically, the system operates as 
follows. Eaeh individual item or component is assigned a unique register 
number which is used for all WIPICS transactions. Each time the status 
or location of the component changes, a transaction updates the computer. 
The system provides immediate audio response (pre-recorded voice) for 
simple questions and printed reports for more detailed questions, A 
query as to the location or status of a component or part is answered 
with current information on that item. 

C, PRIOR THESIS WORK 

The Management Systems Development Office (MSDO) at North Island 
and NARFNI have been charged by higher authority with determining the 
cost-effectiveness of WIPICS as one faetor in justification of continuing 
the use of WIPICS at NARFNI and possible installation of WIPICS at other 
rework facilities. In conjunction with this, a project is continuing at 
the Naval Postgraduate School to develop a methodology for auditing cost- 
effeetivenss analysis of major technological changes. This methodology 
is then to be applied to WIPICS at NARFNI, 

^ Touch Tone (R) is a registered trade mark. 
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Three recent graduates of the Operations Research Curriculum at the 
Naval Postgraduate School each wrote a thesis as part of the above 
project. Spooner [2] wrote on n Evaluation of a Technological Change in 
Production”. He outlined the methodology one might use in determining the 
cost-effectiveness of WIPICS, but he did not actually apply his method- 
ology. In theory, his work takes one from production functions to cost 
functions to product transformation curves. Bradley [l] used actual 
production data taken from ”pre-WIPICS” production records to estimate a 

i 

continuous Cobb-Douglas production function for the aircraft rework, 
engine repair, and component repair programs. Meyers [3] constructed a 
linear economic model for the aircraft rework and engine repair programs 
which predicts resource requirements and costs for several alternative 
objective functions. 

D. SCOPE OF THIS STUDY 

Bradley [l] did not carry his research to the next step, which, 
following Spooner’s methodology, is to estimate a cost function for 
each program of the NARF. There is also evidence which indicates the 
the existence of certain flaws in using Bradley's estimated Cobb-Douglas 
production functions to directly construct cost functions. This is 
cohered in more detail in the next chapter. 

The purpose of this study is to estimate a cost function for the 
aircraft rework and engine repair programs at NARFNI through multiple 
regression analysis of actual production data at NARFNI from March, 1970, 
to August, 1971. This is the second step in Spooner's three step method- 
ology for determining the cost-effectiveness of WIPICS, Using methods of 
reduced form estimation, a cost function will be estimated from a 
Constant Elasticity of Substitution (CES) production function and a 
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Cobb-Douglas production function. The two cost functions will then be 
compared to determine which one more accurately predicts total production 
costs. Finally, a methodology will be outlined for comparing the best 
cost equation and the Linear Programming (LP) Model to the data and to 
each other. 
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II. THE PROBLEM 



A. APPLICATION OF THE COBB-DOUGLAS MODEL 

Spooner [2, p. 33] discusses a computerized prorating program designed 
to prorate the raw data provided by NARFNI for each engine and aircraft 
that underwent repair during a specified period. Assuming seven day work 
weeks and equal distribution of labor hours and material dollars on a 
daily basis, the program prorates NORM, direct labor dollars, direct labor 
hours, direct material costs, and overhead costs for each individual 
engine and aircraft. 

Using aggregated prorated data obtained from the prorate program, 

Bradley [l] estimated Cobb-Douglas production functions for aircraft 

rework and engine repair programs at NARFNI as follows (the equations are 

presented in logarithmic form to facilitate display of standard error of 

2 

each coefficient, shown in parenthesis, and overall R , shown at the 
right of each equation) ; 

Aircraft 

AnAPH = 4.04 + ,38104AnAPL - . 237l4!nAPM + ,88208inN R 2 = .920 

(.044) (.018) (.044) 



Engines 

AnAPH = -1.898 + ,92210£nAPL + . 25135£nAPM R 2 = .981 

(.009) (.011) 

APH is aggregate prorated hours of norm, APL is aggregate prorated hours 
of labor, APM is aggregate prorated material dollars, and N is number of 
jobs in shop. 
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From Naval Air Rework Facility Production and Planning Notices 

dated 17 June 1971, and 26 July 1971, forecast workload requirements for 

engine and aircraft programs during third quarter fiscal 1972 were obtained. 

This information was then used in Bradley* s models to predict quarterly 

direct labor requirements and total cost. One result of comparing actual 

and predicted data appears in Figure 2. Here an isoquant was plotted for 

the aircraft rework program, holding N constant and plotting labor versus 

material. An interesting result immediately became evident - all pro-* 

/ 

duction data points fell in the upper left hand portion of the curve. 
Furthermore, when the isocost line was plotted for the program (using 
prices of labor and material from Meyers* thesis), the theoretical point 
of most efficient operation, which is the point of tangency between the 



LABOR 

Figure 2. Plot of Cobb-Douglas Model (Aircraft) 

isoquant curve and the isocost line (Point A), fell far out of the region 
of actual operation. This implies one of two things; either NARFNI is 
operating very inefficiently, or the Cobb-Douglas production function, 
as Bradley estimated it, is not a good description of the actual process 
by which managers decide on quantities of inputs to use in the repair 
programs at NARFNI. Given the number of variables which enter into the 
production process at NARFNI that are not included in Bradley* s model, 
i.e. managerial decisions and constraints on such items as work flow, 



allocation of resources, and overtime hours permitted, the latter reason 
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r Isocost Line 



given above appears to be more realistic than the former. All of the 
estimation conducted in this study was carried out under the assumption 
that NARFNI was in fact operating efficiently. 

B. NEED FOR A COST MODEL 

The problem is one of estimating a realistic cost function for NARFNI 
as the second step in Spooner T s three step process for determining the 
cost-effectiveness of WIPICS. This cost function can then be used to 
complete the third step, which is formulation of product transformation 
curves . 

Because of the problems encountered in application of Bradley T s 
Cobb-Douglas production functions, it was decided to use methods of 
reduced form estimation to obtain a cost function and provide estimates 
of the coefficients in the CES and Cobb-Douglas models. Another reason 
for not using Bradley 1 s model was the desire to include penalty costs, 
as determined by Meyers [3], in a cost function as a measure of pipeline 
cost per unit. Penalty costs are the dail}' costs to the Navy for having 
an aircraft or engine in the rework cycle. Penalty costs for aircraft 
are computed by dividing procurement cost by expected life. Penalty 
costs for engines are determined by dividing procurement cost by ten 
years. By including an index of the price of labor (prorated labor 
costs divided by prorated labor hours) and an index of pipeline costs 
(aggregated penalty costs divided by daily number of items in shop), 
predicted costs will be responsive to increases in wages as well as 
increased sophistication of the items undergoing rework. 
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III. DERIVING COST FUNCTIONS 



A. THEORY 

1 . General 

The theory behind the derivation of cost functions in this 
chapter is based on three very important assumptions. First, we assume 
that it is the objective of the NARF to minimize costs. The second 
assumption is that the NARF has no control over prices of labor or 
penalty costs per item per day. Third, we assume that output is fixed 
at any desired level and the NARF does an imperfect job of minimizing costs 
for a given level of output. The reason for this is that when the NARF 
makes a decision it doesn't know exactly what the prices will be, and it 
does not know until an item arrives exactly how extensive a rework will 
be required. According to Nerlove [4, p. 107], the fundamental duality 
between costs and production functions assures us that the relation 
between the cost function, obtained emperically, and the underlying 
production function is unique. 

2 . CES Production Functio n 

The general form of the CES production function is as follows 
[5, p. 284]: 

-o -d ~ a/p 

Y = 7(6X 1 P + (1-6)X 2 P ) 

where y is an efficiency parameter, p is a substitution parameter, and 
6 is a distribution parameter, a is a measure of returns to scale, i.e., 
if o = 1 , there are constant returns to scale, and if a > 1 or a < 1 
there are increasing or decreasing returns to scale respectively. Y is 
production output, while X^ and are inputs to production. 
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What follows is the derivation of a cost function from the production 
function and marginal productivity relations. See Walters [5, p. 284] 
for a more detailed explanation. 



min C = P^X^ + 



( 1 ) 



s.t.Y = 7(6X 1 “ P + (1 - 6)X 2 ' P )" a/p 



where C is production cost, P is the price of input X- , and P is the 

i L £• 

price of input X 2 « Forming the Lagrange Equation and taking partial 
derivatives , 



Z = P 1 X 1 + P 9 X 2 - X[Y - 7(6X~ P + (1 - 6)X 2 “ P )" ct/p ] 



( 2 ) 



H~ = ~ 7(8X 1 ' P + (1 - 5)X 2 “ p ) ( ' Cj/p) " 1 (l - 6)(-p)X 1 " P ” 1 = XP 1 (3) 



_ 22. 
BX 0 " p 



-p. (-a/p) -1 



, -P-1 = 



7(5X 1 K + (1-6)X 2 (5) (-p)X 2 ^ = XP 2 



(4) 



Forming the ratio of equations (3) and (4) , 



(1 - 6)X 1 ' P " 1 



6X. 



-p-i 



(5) 



Thus , 



X, 



6P, 



(1 - 6)P 2 



-P-1 



X 0 



( 6 ) 



where |i is a random error term due to imperfect minimization of costs. 
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Now, 



-P-1 



X„ 



(1 - §)P 2 



(7) 



Taking natural logarithms of both sides of (7) 



In (X l /X 2 ) = a + b^n(P 2 /P 1 ) + In p, 



( 8 ) 



where 



-P-1 



1-6 



and b= 



•P-1 



Substituting back into the cost equation. 



X n = 



2 P 1 (X 1 /X 2 ) + P 2 



(9) 



The production function is now 



Y = y([6 + (1 - 6) (x ] /x 2 )~ p ]x 2 )" cr/p 



( 10 ) 



Solving for X,. 



(r 



yp/a 



6 + (1 - 6)(X 1 /X 2 ) 



, x“P 



1/-P 



( 11 ) 



Taking natural logs of both sides of equation (11), 



AnX„ 



a 



In 



y 



+ u M+ (i 



6 + (1 - 6)(X 1 /X 2 ) 



( 12 ) 
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To simplify handling equation (12), let 



XnX 2 = |>] + 3XnY + (u>)[d] 



( 13 ) 



The final form of the cost equation is 



C 




(14) 



where e is a random error term. 



3 . Cobb-Douglas Production Function 



The general form of the Cobb-Douglas production function is 



[5, p. 275]: 




where Y is production output, and are production inputs, and O' and 
(3 are elasticities of production with respect to X^ and X 2 respectively. 



Once again the specification of the production function is changed so that 
Y is non-random, but X^ and X 2 sre random variables, due to the fact 
that the NARFNI does an imperfect job of minimizing costs for a given 
level of output and has no control over prices. With this in mind, a 
cost function is derived as follows: 



min C = P ][ X 1 + P 2 X 2 



(15) 




Forming the Lagrangian Equation and taking partial derivatives, 




( 16 ) 
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H: ■ p i • toAx rv “i ■ ° 



(17) 



Hr “ p 2 - “ 2 ■ 0 



(18) 



|f - 0 



(19) 



where the p ! s are random variables associated with error (again, due to 
imperfect minimization of costs). Taking the ratio of equations (17) 
and (18), 



c,'X 2 



li 



( 20 ) 



Solving for inputs X.^ and X 2 , the following cost function is obtained: 



C = 



Y a 
A [ B 



+ 1 



e. + 

O' 



va a p 
P 1 P 2 



- 1 1/(Q'+B) 



f‘ 



-_!L-/gli + A (21) 

cy + B l a J 



B, PRODUCT TRANSFORMATION CURVES 

Spooner [2, Appendix E] contains a complete theoretical derivation 
of product transformation curves from cost functions. Basically, a 
product transformation curve is obtained by rewriting a cost function in 
terms of output. The result of the substitution is a nonlinear rela- 
tionship between two or more outputs. By looking at two outputs at a 
time, setting the others equal to a constant, and setting 'cost equal to 



20 



a constant, graphical analysis can be applied to the overall cost function. 
After plotting one output against the other, the shape of the product 
transformation curve will indicate increasing, constant, or decreasing 
returns to scale (see Figure 3) . 



Output 



Y 



2 




Figure 3, Product Transformation Curves, Returns to Scale 
Product transformation curves may also be used as one way of compar- 
ing the effects of WIPICS on production at NARFNI. A comparison could 
be made between before WIPICS and after WIPICS transformation curves, 
keeping the total budget constant and keeping prices of resources and 
inputs the same for the two periods. One should be very careful in 
evaluating the results of such a comparison because of the many factors 
which enter into production at NARFNI but do not appear in a cost function, 
particularly management decisions. 
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IV. ESTIMATING COST FUNCTIONS 



A. PRORATING WITH PENALTY COSTS 

The raw data on each aircraft and engine job completed during the 
period February, 1970 to August, 1971, which appears in Meyers' thesis 
[3, Appendix A] was run through the computer prorate program in Spooner's 
thesis [2, Appendix B] in order to provide daily aggregated, prorated 
data which could be used in a linear regression. In addition, penalty 
costs for each type of engine and aircraft (see Table I and Table II) 
were incorporated into Spooner's prorate computer program so that these 
costs could be used in estimating a cost function for NARFNI, The codes 
shown in Tables I and II are as defined by Meyers [3, p, 16], The method 
used to incorporate these costs into the program was simply to read the 
aircraft or engine code from the raw data card for all items that were in 
shop each day, and then the appropriate daily penalty costs were aggregated 
for each day of the period under observation. The aggregated prorated 
data used to estimate the cost functions in this thesis appear in 
Appendix A (aircraft) and Appendix B (engines). Spooner's prorate 
programs, rewritten to include penalty costs, appear in Appendix D 
(aircraft) and Appendix E (engine) . 

The interest in penalty costs stems from the fact that this appeared 
to be the only way to quantify unit pipeline costs, and it was considered 
desirable to include pipeline costs in any cost model of NARFNI . 
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TABLE I 



DAILY ENGINE PENALTY COSTS 



CODE 


TYPE ENGINE 


PENALTY COST 


51 


T58-GE-8B 


17.81 


53 


T58-GE-8B 


17.81 


55 


T58-GE-8B 


17.81 


57 


T58-GE-8B 


17.81 


81 


T58-GE-8B 


17.81 


82 


T58-GE--8B 


17.81 


83 


T58-GE-8B 


17.81 


84 


T58-GE-8B 


17.81 


85 


T58-GE-8F 


17.81 


86 


T58-GE-8F 


17.81 


87 


T58-GE-8F 


17.81 


89 


T58-GE-8F 


17.81 


63 


T64-GE-6B 


34.25 


65 


T64-GE-6B 


34.25 


66 


T64-GE-6B 


34.25 


67 


T64-GE-6B 


34.25 


69 


T64-GE-413 


45.25 


56 


T56 ALL MODELS 


27.01 


71 


J57 ALL MODELS 


57.53 


72 


J57 ALL MODELS 


57.53 


73 


J57 ALL MODELS 


57.53 


74 


J57 ALL MODELS 


57.53 


75 


J57 ALL MODELS 


57.53 


91 


J79-GE-8B 


61.64 


92 


J79-GE--8C 


61.64 


95 


J79-GE-10 


49.39 
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TABLE II DAILY AIRCRAFT PENALTY COSTS 



CODE 


TYPE AIRCRAFT 


PENALTY COST 


10 


C-2A 


280.29 


11 


E-2B 


280.29 


21 


F-4J 


85.34 


22 


F-4B 


94.38 


23 


RF-4B 


77.02 


25 


F-8J 


66.05 


26 


F-81I 


68.59 


27 


RF-8G 


47 . 48 


31 


CH-3B 


43.76 


32 


RH-3A 


51.64 


33 


SH-3A 


43.76 


34 


SH-3A/G 


36.43 


35 


SH-3D 


36.43 


41 


CH-46A 


37.08 


42 


CH-46D 


37.08 


43 


CH-46F 


37.08 


44 


UH-46A 


28.56 


45 


UH-46D 


28.56 


48 


CH-53A 


41.46 


49 


CH-53D 


41.46 
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B. REGRESSION PROCEEDURES AND RESULTS 
1. CES Proceedures 

. For the CES production function a cost function was estimated as 
follows (all values used were taken from the prorated data, cut on each 
end to eliminate start up and shut down effects as per Spooner [2, p fl 33]) 

min C = PI (I) + PL (L) 

-a/p 

s.t.N = y (SL~ P + (1 - 6)I~ P ) 



where 



PI 



penal t y co s t 
# of items in shop 



PL 



direct labor cost 
direct man hours 



I - number of items in shop 



L = direct man hours 



N = NORM 



Following the proceedure as set forth in the previous chapter, 
equations (1) through (8), the first stepwise linear regression was run 
on equation (8) of Chapter 3. This equation is listed below, with in \h 
representing the natural log of the error term: 



in(l/L) = a + bXn (PL/PI) + in pi 

1 



where a 



.(-JL.P 
V 1 - v 



and b = 



(i) 



The error was assumed to be 



P + 1 

normally distributed with a mean of zero and constant variance. 
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The regression on equation (1) provided an estimate of (i/L), which 

A 

will be indicated as (i/L), for use in the following equation: 



L* 



C 

PI (I/L) + PL 



(2) 



which corresponds to equation (9) in the previous Chapter. Here C equals 
direct labor cost plus overhead cost plus penalty cost. The regression 
on equation (1) also provided the parameters for obtaining the quantity 

A ' 

D in equation (4) below. The final equation on which a second stepwise 

\ 

linear regression was run was (see equation (12), Chapter III): 
jJnL* = 

where e is a random error term. Simplifying the quantities in equation 
(3), let 

A 

j&nL* = a + b AnN + c[D] + e (4) 

The second regression provided estimates of the coefficients of the 
independent variables. 

All regressions in this thesis were run using the SNAP/IEDA statistical 
package available at the Naval Postgraduate School Computer Center. 

2 . CES Regression Pvesults 

The results of the first CES regression for the aircraft and engine 

2 

programs are presented in Tables III and IV respectively. The low R 
(.05) in the aircraft regression is an indication of possible problems 

in applying the CES production function to the production process at 

2 

NARFNI (R is the square of the multiple correlation between the dependent 
variable and those independent variables which were included in the 



— in y 

a f 



+ (J)^N + (~) 



* “P 

6 +(1 - 6) (I/L) 



+ e 



( 3 ) 
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TABLE III. FIRST CES STEPWISE REGRESSION OF AIRCRAFT DATA 



STEP 



Entered 



Previously Entered 



R 



Standard Error of Dependent Var. 



1 

PL/PI 

.054 

.036 



Coefficients 
Constant 
PI /PL 



- 4.010 

.233 



Standard Error of Coefficients 
PI /PL 



.058 



CV 



■.0078 
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TABLE IV. FIRST CES STEPWISE REGRESSION OF ENGINE DATA 



STEP. 1 

Entered PL/PI 

Previously entered 

R^ .440 

Standard error of dependent var, .074 

Coefficients 

Constant -0.219 

PI/PL 1.324 

Standard Error of Coefficients 

Pi/ PL .079 

CV - .0266 
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regression at that step -- it measures the proportion of variation in 

the dependent variable that is explained by the independent variable) . 

2 

Rao and Miller [6, p. 16] point out that a high R may imply the 

appropriateness of a regression equation for explaining the movements 

2 

of a dependent variable, but a low R does not necessairly imply that the 

regression equation is inappropriate . Another possible problem 

appeared when the value of 6 was computed. It is specified by Walters 

[5, p. 284] that 0 ^ 6 £ 1 in a CES production function. This means that 

the expected value of the constant in the first regression runs was a 

positive number less than one. This was not the case for either progiavn. 

As a matter of interest, the coefficient of variation (CV) , which is 

standard error divided by mean of the dependent variable, is listed in 

all tabled regression results in this Chapter, It is a measure of 

dispersion of the dependent variable in comparison to the mean due to 

lack of a perfect fit with the data. 

Looking at the results of the second regression of the aircraft data 

2 

(Table V), the extremely high R (,988) indicates a good fit with the 

data. Both independent variables, N and D, were found th have significant 

coefficients after examination of their respective standard errors (when 

computation of the standard deviation is based on the estimate of the 

variance rather than the variance itself, it is called standard error). 

In the engine program (Table VI) , both independent variables were 

again found to have coefficients significantly different from zero. 

in this case (.798) was not nearly as high as it was for the aircraft 

program, but it is certainly acceptable. 

The value of p, computed from the coefficient of D in each case 

(the coefficient of D equals - ~) , shows p to fall well within the 

p 
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TABLE V SECOND CES STEPWISE REGRESSION OF AIRCRAFT DATA 



STEP 

Entered 

Previously entered 

2 

R 

Standard error of dependent variable 



1 

N 



.977 

.012 



2 

D 

N 

.988 

.009 



Coefficients 

Constant 

N 

D 



1.596 2.063 

.897 .880 



- 2.388 



Standard Error of Coefficients 

N 

D 

CV 



.008 .006 
.149 

.0006 



30 



TABLE VI SECOND CES REGRESSION OF ENGINE DATA 



STEP 1 2 

Entered N D 

Previously entered N 

R 2 .775 .798 

Standard error of dependent variable .066 .062 

Coefficients 

Constant .753 2.344 

N .983 .929 

D 4.566 

Standard Error of Coefficients 

N .028 .028 

D .721 

CV .0075 
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expected region. This value was ,419 for the aircraft program and -,219 
for the engine program. The expected value of p is a number between minus 
one and infinity. As p approaches minus one the result is a flat isoquant, 
and as p increases in a positive direction the elasticity of substitution 
becomes smaller and smaller (essentially this is the case of fixed 
coefficients). For a more detailed explanation of this, see Walters 
[5, p. 286]. 



The final estimated cost functions were: 



/S. 



(aircraft) 



(engines) 



C- (PL + PI(I/L»e 2 - 06 V 880 e- 2 ’ 388 < D ) 
C- (PL + PI(I?L» e 2 - 34 V 929 e 4 ' 566 ^ 



( 5 ) 

( 6 ) 



3. Cobb-Douglas Proceed ures 

For the Cobb-Douglas production function, a cost function was 
estimated using the general theory presented in equations (15) through 
(21) in Chapter III, The specific equations used in the estimation 
proceedure were: 



min C = PL(L) + M + PI(I) 
s.t.N = AL Q I^M y 



All variables are as previously defined, and M is material costs (again, 
prorated data was used in all cases) . A stepwise linear regression was 
run on the following equation: 



log C = a 



Ol + 3 + y 



log N 4* 



a 



H- ~ 



a + (3 +y 

p 

a? + (3 + y 



log PL 



log PI + e 



(7) 
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where a 



log A 




o' + P + y 



1 




Again, the error term is assumed to be normally distributed with a mean 
of zero and constant variance. The regression provided estimates of the 
coefficients of N, PL, and PI, 

4 , Cobb -Doug las Reg r ession Results 



in Tables VIT and VIII for the aircraft and engines programs, respectively. 

Examination of this tabulated data prvides a good indication that the 

cost functions estimated from the Cobb~Douglas production function may 

turn out to be valid and accurate cost functions „ 

In the aircraft data regression, all coefficients were significant, 

2 

and the value of R (.960) was an indication of a good fit with the 

data. The regression on the engine data showed only NORM to have a 

significant coefficient. The coefficient of PL(~„027) had a standard 

error of .113, and the coefficient of PI(.101) had a standard error of 

.051, indicating that neither coefficient was significantly different 
2 

from zero. R for the engine data was ,848. 

The resulting cost functions for the aircraft and engine programs 
were : 



The results of the regressions run on. equation (!) are presented 



(aircraft) C = 10 



• 42 V 751 PI. 2 - 629 PX " 275 



( 8 ) 



(engines) C 



lo 1 * 749 ^ 860 



(9) 
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TABLE VII COBB-DOUGLAS STEPWISE REGRESSION OF AIRCRAFT DATA 



STEP 


i 


2 


3 


Entered . 


N 


PL 


PI 


Previous ly entered 




N 


N 








PL 


R 2 


. 645 


.949 


.960 


Standard error of dep. var. 


.024 


.009 


.008 


Coefficients 








Constant 


1.678 


-0.049 


0.421 


N 


0.824 


0.769 


0.751 


PL 




2.453 


2.629 


PI 






-0.275 


Standard error of Coefficients 








N 


.036 


.014 


.012 


PL 




.060 


.056 


PI 






.030 


CV 


.0049 


.0018 


.0016 
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TABLE VIII COBB-DOUGLAS STEPWISE REGRESSION OF ENGINE DATA 



STEP 


i 


2 


3 


Entered 


N 


PI 


PL 


Previously entered 




N 


PI 








N 


2 

R 


.848 


.849 


.849 


Standard 6rror of dep c var. 


.020 * 


.020 


.020 


Coefficients 








Constant 


1.749 


1,539 


1.554 


N 


.860 


.876 


.876 


PI 




.097 


.101 


PL 






-.027 


Standard error of coefficents 








N 


.019 


,021 


.021 


PI 




.047 


.051 


PL 






.113 


CV 


.0044 


.0044 


.0044 
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V. VALIDATION 



A. APPLICATION OF THE COST MODELS 

The prorated data was broken down into nine consecutive 30 day 
intervals beginning with julian date 0191 and ending with julian date 
1095. Average daily values of NORM, total costs (labor plus material 
plus penalty), PL, and PI were computed for each period. This was done 
for both the aircraft and engine programs (see Appendix C) . 

For validation of the CES cost functions, the values of NORM, PL, 

A ^ 

PI, i/L, and D were computed for the nine 30 day periods and then used 
in equations (5) and (6) in Chapter IV. The results of these runs, 
shown in Table IX, are evidence that the CES cost function can in fact 
be used to predict costs at NARFNI (remembering that actual costs shown 
here are average daily prorated direct labor cost plus overhead cost 
plus penalty cost). The worst prediction occurs in the first engine 
prediction, where the difference between predicted and actual cost is 
only approximately four thousand dollars out of thirty eight thousand, 
or less than 107 o . The percent of error for each program as shown in 
Tables IX and X was computed by subtracting actual cost from predicted 
cost and dividing by actual cost. The same independent variables used 

A A 

above, less i/L and D, were then fed into the estimated Cobb-Douglas 
cost functions which appear at the end of Chapter IV (equations (8) and 
(9)). Actual and predicted costs for the aircraft and engine programs 
are listed in Table X for each of the nine 30 day periods. An examination 
of the data in Table X shows the Cobb-Douglas cost functions to also 
be predicting costs with fairly good accuracy (again, a reminder that 
actual costs here are the sum of direct labor, material, and penalty costs). 



TABLE IX ACTUAL AND PREDICTED TOTAL COSTS (CES) 





AIRCRAFT 






PREDICTED 


ACTUAL 


% 


ERROR 


96,896.19 


98,561.24 


- 


1.69 


116,596.25 


117,882.93 


“ 


1.09 


113,573.44 


112,113.25 




1.30 


116,616.94 


115,304.52 


+ 


1.14 


119,525.25 


118,598.68 




.78 


110,645.50 


109,015.75 


+ 


1.50 


113,536.06 


113,658.38 


- 


.08 


121,063.69 


122,274.12 




.99 


114,884,63 


115,261.18 


- 


.33 




ENGINES 






PREDICITED 


ACTUAL 


7o 


ERROR 


34,439.70 


38,230.57 


- 


9.93 


36,684.21 


39,975.95 


- 


8.25 


30,655.29 


32,061.61 


« 


4.38 


36,398.64 


36,459.35 


- 


.17 


36 , 644.64 


36,471.29 


+ 


.48 


29,525.21 


30,294.42 


- 


2.55 


35,207.98 


34,942.32 




.76 


33,748.04 


32,864.77 


+ 


2.70 


29,846.05 


29,806.54 


4- 


.13 
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TABLE X ACTUAL AND PREDICTED TOTAL COSTS (COCB-DOUGLAS) 





AIRCRAFT 






PREDICTED 


ACTUAL 


% 


ERROR 


67,304.50 


68,452.51 


- 


1.68 


79,858.06 


80,122.31 


- 


.33 


77,338.25 


74,926.64 




3.26 


78,968.56 


78,632.10 


+ 


.43 


81,799.38 ' 


82,816.17 


- 


1.23 


77,561.00 


77,510.11 


+ 


.06 


82,942.81 


82,667.20 


+ 


.33 


89,536.81 


89,337.36 


+ 


.22 


84,495.63 


85,059.84 


- 


.67 




ENGINES 






PREDICTED 


ACTUAL 


% 


ERROR 


41,720,43 


45,946.33 


- 


9.20 


44,015.41 


45,183.24 




2.59 


37,536.84 


36,286,71 


+ 


3.45 


43,974.30 


41,814.51 


+ 


5.16 


43,680,50 


41,564.47 


4* 


5.09 


34,941.77 


33,638.37 


+ 


3.90 


41,488.91 


41.437.46 


+ 


.12 


39,378.77 


39.130.91 


+ 


.63 


35,189.98 


36,008.00 


- 


2.27 
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B. PREDICTION INTERVALS 



A prediction interval can be obtained for the Cobb-Douglas aircraft 
cost function and the engine cost function using methods -outlined in 
Theil [7, p. 135], This method proceeds as follows. First, return the 
cost equation to logarithmic form. Then let, 

Y* = W ' C + e* 



where Y* is a single observation of cost in log form, W* is a one row 
vector, with each column containing given values of the independent 
variables, C is a one column vector with each row being a coefficient of 
one independent variable, and e* is the random error term with a mean of 
zero. Let Z be a matrix of observations on which the regression was 
run. The first column is all ones, and the remaining columns are data 
columns. There are as many rows as there are numbers of observations. 



i • e * > 



1 



N, 



PL, 



PI 



1 



Z = 



N 



PL. 



N 



PI. 



W * C - Y* 

Now =y=— has a student's t distribution with n-k degrees 

VS (1 H- W' (Z'Z)”^) 

of freedom. Here n equals number of observations and k equals number of 

2 

independent variables in the cost equation, S is an estimate of 
variance. Selecting a 1 - ot confidence interval. 



T ^~ t o'/2 £ t (n-k) ^ L/2 ^ 



1 -<v 



Let 



S 



1 + W'(Z'Z) -1 W , so that 



V = t 



o72 



p[w'c -VS Y* < W'C + v] = 1 - or 



Taking anti logs. 



[10 (H ’ C - V > S C S 10 (H ' C+ v) 



] 



is a 1 - O' prediction, interval for cost. 

C. METHODOLOGY FOR COMPARISON WITH LP MODEL 

A good way to determine the validity of an estimated cost function 
would be to compare actual costs with estimated costs obtained from the 
cost function and estimated costs obtained from Meyers [3] LP model. 

To make the comparison realistic, several important factors must be 
considered. First, the objective function of the LP model must be to 
minimize costs. Second, the values of NORM and PL, as well as penalty 
costs and material costs used in the LP model, must also be used in the 
cost function. In the LP model, the actual values of NORM for the period 
of time under observation should be used in the production vector Y. 

This constrains the model to minimize costs subject to actual values of 
NORM. The activity vector Z need not be changed, but prices of labor 
and material used in the resource vector R must coincide with the prices 
used in the cost function. There remains the task of breaking the data 
down by type of aircraft or engine and the type of work done on that item 
for use in the activity vector Z. 

Ideally, the two methods of cost prediction should each be run over 
several periods of time, either 30 day intervals, or monthly, or perhaps 
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quarterly. This provides enough data points to carry out a meaningful 
comparison of the two models. 

The criteria for determining which cost model is better at predicting 
costs, as outlined below, is sketchy at best. The lack of independence 
between the time periods eliminates most conventional statistical tests 
which might apply. Naturally, if one model is consistently and obviously 
better than the other, no test need be applied. But if one model is 
better than the other at times, and then at times worse then the other, 

t 

the criteria outlined here may aid in determining a "best" model. 

The criteria is simply this; take the difference between actual and 
predicted costs for each time period and square this number. Then sum 
the squared numbers over the periods of time. Doing this for each model 
will provide a "mean squared error 11 (MSE) term for the models, 



MSE 



i = l 



(C 



Pi 



C Ai> 



2 



where 



n = number of periods 



C . = predicted cost for period i 

pi 



C Ai = actual cost for period i 

The model with the smallest MSE should be considered the best cost model. 
Naturally, if the two values of MSE are very nearly the same, there 
remains a problem of determining whether or not they are significantly 
different. This problem will not be addressed in this thesis. 

As a matter of interest, the values of MSE for the CES and Cobb- Douglas 
cost functions estimated in Chapter TV of this thesis were as follows: 
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CES aircraft 



1,486,812 
CES engines 3 , 180,7 10 
Cobb-Douglas aircraft 1 , 083 , 037 
Cobb-Douglas engines 3 ,584 , 340 
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VI , CONCLUSIONS AND FUTURE STUDY 



A, CONCLUSIONS 

The cost functions, estimated from the CES production function were: 
(aircraft) C = (PL + PI(I/L)) e 2 *° 63 N* 880 e " 2 * 388 (D) 

(engines) C = (PL + PKI/L) ) e 2 ’ 3 %' 92 V' 566 ™ 



These cost functions were concluded to be a valid means of predicting costs 
at NARFNI after several validation runs showed predicted and actual 
average prorated daily costs to compare favorably. 

The cost functions estimated from the Cobb-Douglas production 
function were : 



(aircraft) C - 10' 421 N' 751 PL 2 ' 629 Pl'- 275 



(engines) 



c - lO 1 - 749 *- 860 



The coefficients of PL and PI were not significant in the engine program. 
A possible explanation for the relatively high coefficient of PL and 
the negative coefficient of PI in the aircraft program is that as the 
cost per unit of items being repaired increases, penalty costs increase, 
and due to the complexity of the more expensive items, number in shop 
decreases. Also due to the complexity of more expensive items being 
repaired, the cost of labor increases rapidly as more skilled tradesman 
are required. 

The cost functions derived from the Cobb-Douglas production function 
were also concluded to be a valid means of predicting costs at NARFNI. 
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Predicted and actual average daily prorated costs compared favorably when 
examined over nine 30 day intervals. 

Based on the methodology outlined in this thesis for comparing cost 
prediction models, the cost functions estimated in this thesis predict 
costs with approximately the same accuracy. 

B. FUTURE STUDY 

A prime area for future study would be to compare the LP model, the 
Cobb-Dougias cost equation estimated in this thesis, and actual total 
production costs at NARFNI for several selected periods of time, using 
the methodology outlined in Chapter IV of this thesis. This would 
provide a more firm idea of just how accurate and useful the Cobb-Douglas 
cost function really is. 

Another area that could be improved is the method used to determine 
penalty costs. Rather than simply dividing procurement cost by an expected 
life for each type aircraft or engine, a more realistic method to use 
could be as follows: 



Penalty Cost 



Cost of all a/c, allowing Cost of all a/c, without 

rotable spares rotable spares 



In other words, the true cost to the Navy for having an aircraft at 
NARFNI is the cost of an aircraft that had to go to the fleet to replace 
it (rotable spares) , Dividing this figure by the average length of time 
an aircraft or engine is undergoing rework would provide a daily penalty 
cost . 

Finally, some method of including management decisions and constraints in 
a production function must be devised. New equipment, improved proce- 
edures, increased worker efficiency as a result of labor force cuts. 
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hiring and firing policies, and constraints on such items as overtime 
hours permitted, all enter into the production process at NARFNI, but 
finding a method to quantify these areas and incorporate them into 
a production function is a most difficult problem. 
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SUMMARY OF DAILY STATISTICS OF NARr (ENGINES) 
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APPENDIX C 



AVERAGE DAILY AGGREGATED DATA FOR NINE 30 DAY PERIODS 



NORM 



AIRCRAFT 

7 143 „ 99 
8735.46 
8400.21 
8602.12 
8746.54 
7874.49 
8071.96 
8637.05 
8160.86 



ENGINES 

2181.62 

2321.78 

1929.38 

2319.26 

2301.25 

1775.17 

2167.55 

2039.90 

1789.84 



PL 



AIRCRAFT 

5.99 
6 , 05 
6.07 
6.09 
6.15 
6.23 
6.31 
6.36 
6.33 



ENGINES 

5.92 

5.97 

5.87 

5.87 

6.02 

6.15 

6.15 

6.25 

6.13 



PI 



AIRCRAFT 

86.75 

87.33 

91.01 

92.89 

93.93 

96.81 

91.69 

90.06 

91.26 



ENGINES 

40.58 

40.50 

41.07 
41.17 
40.68 

42.08 
40.76 
40.97 
42.84 
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AIRCRAFT PRORATE PROGRAM 
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